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Although of primary importance in the mechanistic understanding of Zircaloy hydriding, practi- 
cally nothing is known on the transport properties of hydrogen in zirconia. In this frame the muon, 
which can be considered as a light hydrogen nuclide, can be used as an analogue and its behavior in 
zirconia and Zircaloy corrosion layer may provide more insight to understand the behavior of hydro- 
gen in these phases. Preliminary muon spin relaxation (/xSR) measurements on several monoclinic 
zirconia samples, including a Zircaloy corrosion layer, have been performed. From the observed 
muon depolarization rate, the muon diffusivity in bulk monoclinic zirconia can be extracted and is 
found comparable to that of recently reported proton diffusivity. 
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INTRODUCTION 

The study on the behavior of hydrogen in zirconia for 
understanding its migration mechanism in the Zircaloy 
corrosion layers is justified as stated in specific studies 
and as reported recently (see Ref. 0): 

Virtually, nothing is known about the form of 
hydrogen in Zr02 or how it migrates through 
it. (...) The arguments have always been over 
precisely what changes are occurring in the 
oxide during this incubation time, and how 
the hydrogen migrates through the oxide. 

An atomistic picture of interstitial hydrogen in zirconia 
involves, as it does for hydrogen in dielectrics, determina- 
tion of lattice sites, configurations adopted and the mech- 
anisms of diffusion. An additional dimension is added in 
the study of muons in materials, namely the possibility 
of examining different states of the defects centers - posi- 
tive, neutral, negative - respectively for the diamagnetic 
proton, the paramagnetic hydrogen atom and diamag- 
netic hydride ion. Muon spectroscopy has been success- 
ful in characterizing all these states, determining their 
local structures, their different stabilities, mobilities and 
interactions with charge carriers 0|. In addition charac- 
terization of trapping centers such as vacancies may be 
performed. Such similar studies could therefore be em- 
ployed to get more insight on the behavior of hydrogen 
in for zirconia, as the muonium {Mu = /i+ + e~) can be 
regarded as an isotopic analogue of the hydrogen atom, 
with a positive muon replacing the proton. 

In metals, the implanted muon comes at rest at an 
intersticial position and is in a diamagnetic state. Af- 
ter implantation, it may diffuse between interstitial sites 
and/or towards trapping spaces (defects, vacancies, for- 
eign atoms). By /iSR technique, information about the 
local magnetic field distribution at the muon site can be 
obtained. Such distribution will be sensitive to the dy- 
namics of the local moments and/or of the muon itself. 



Consequently muon diffusion in metals may be investi- 
gated and numerous studies are reported for elemental 
metals (see for example Ref. 3]). Alternatively, in di- 
electrics, part of the implanted muons can capture an 
electron and form a muonium pseudo-atom. The Muo- 
nium properties, including the nature of trapping sites 
and its diffusion, can in turn be studied by /iSR, as suc- 
cesfully demonstrated in a variety of studies (for a com- 
prehensive review on muon and muonium diffusion see 
Ref. H). 

As a complement of recent studies using x-ray absorp- 
tion 01 to gain more insight on the influence of oxy- 
gen deficiency on electronic and local structure of mono- 
clinic zirconia, the present /iSR investigation aims to ob- 
tain complementary information on the muon diffusion 
on monoclinic zirconia (powder and polycrystalline sam- 
ple). Moreover, preliminary data on a Zircaloy corrosion 
layer are also reported. 

EXPERIMENT 
Muon Spectroscopy 

The muons used during /zSR studies are basically 100% 
polarized and are implanted one at a time into the sam- 
ple, where they thermalize within few picoseconds. Dur- 
ing his lifetime (ca. 2.2 /zs), the muon can diffuse in the 
zirconia material as free (diamagnetic) muon or possibly 
as muonium. When the muon undergoes its weak decay 
a positron, among others, is produced which is prefer- 
entially emitted along the direction of the muon spin at 
decay time. For a /iSR experiment, the direction along 
which the decay positron is emitted and the elapsed time 
interval between the muon implantation and the decay 
are determined. After ~ 10^ muons have been observed, 
the obtained time histogram of the collected intervals has 
the form: 

Nit) = No exp(-i/T^)[l + AG{t)] + Bg , (1) 
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where Bg is a time-independent background, Nq is a nor- 
malization constant, is the muon Ufetime and the ex- 
ponential accounts for the muon decay. AG{t) is often 
called the /iSR signal. A is the asymmetry parameter, 
which will depend on the experimental setup. G(t) re- 
flects the normalized muon-spin S autocorrelation func- 
tion G{t) = (S(i) • S(0))/S'(0)^ which depends on the 
average value, distribution, and time evolution of the 
internal fields sensed by the muons, and therefore con- 
tains all the physics of the magnetic interactions of the 
muon inside the sample Q. The present /iSR measure- 
ments were obtained with the so-called transverse-field 
technique where a weak magnetic field is applied trans- 
verse to the initial spin direction. The form and width of 
the field distribution sensed by the muon is therefore de- 
tected by monitoring the decay of the muon polarization 
(i.e. reflected by the envelope g{t) of the G{t) function). 

The /xSR experiments were performed at the Swiss 
Muon Source (S/iS) of the Paul Scherrer Institute (PSI), 
Villigen (Switzerland). The GPS instrument located 
on the surface muon 7rM3.2 beam-line was used. The typ- 
ical range in matter of the muons available on this beam- 
line is about 150 mgcm~^. Samples were cryo-cooled 
in liquid helium and /iSR was performed as a function of 
temperature between 5 and 300 K. Due to the small depo- 
larization rates observed, the so-called Muon-On-Request 
(MORE) setup was used. This setup dramatically re- 
duces the value of Bg and allows one to extend the usable 
/iSR time window up to 20 /is and therefore detect depo- 
larization rates as low as 0.001 /iS. 



Samples 

Different monoclinic zirconia samples were investi- 
gated: i) a very pure powder sample; ii) a solid pel- 
let sample; iii) a sample of natural monoclinic zirconia 
formed by a mosaic of single crystals. The zirconia pow- 
der was produced by Fluka (pro-analysis) and was com- 
posed by particles in the size range of 1 to 10 fi. The 
zirconia solid pellet was a sample of nanoparticles of 
pure zirconia pelletised and sintered at low temperature 
(kindly offered by the PennState University) as reported 
earlier by Raghavan et al. . Natural zirconia was part 
of a brown rolled pebble comprising baddeleyite (Zr02) 
and zircon from PoQOS de Caldas, Minas Gerais (Brazil). 
The baddeleyite part was 22.5 x 22 mm^ with a thick- 
ness of 1.5 to 3 mm. This sample was kindly provided by 
the Mineralogical Department of the Geological Museum 
of the University of Lausanne (Switzerland). The main 
impurities were Fe and Hf both at the percent level. In 
addition, a Zircaloy (Zry4) sample was also investigated. 
This sample (40 x 10 mm^), corroded in autoclave, was 
covered by a layer of Zr02. This monoclinic zirconia cor- 
rosion layer had a thickness of 100 /im as determined by 
gravimetry (1692 mgdm~^). All samples were tested by 
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FIG. 1: Upper panel: Calculated penetration depth for lep- 
tons and hadrons with a 4 MeV incident energy (see text). 
Lower panel: Comparison between the calculated penetration 
depths of /i"*" (squares), H"*" (triangles) and ^H^ (circles). 



x-ray diffraction and were found to assume the mono- 
clinic structure of zirconia. 

During /iSR measurements, the samples were mounted 
either on special containers (powder sample) or on fork- 
like holders. In both cases, the fraction of muons missing 
the sample and creating a time-dependent background 
signal can safely be neglected. 
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FIG. 2: TF nSK spectra, recorded in an external field of 7 G 
over the full available time window (8 fis - top) and over the 
first 0.5 us (bottom), in monoclinic zirconia at low temper- 
ature (4 K). The signal surviving at long time is associated 
to diamagnetic muons, whereas the rapidly damped signal is 
ascribed to a muonium fraction (see text). 



RESULTS AND DISCUSSION 
Modelling the Muon Implantation 

Penetration and implantation of leptons "a:"'" and 
hadrons "X"*" in zirconia (density 5.9 gcm~'^) were mod- 
eled using the software SRIM2000 (Ref. 11]). First cal- 
culations were carried out for energy of 4 MeV corre- 
sponding to the muon energy available for the real mea- 
surements. The average penetration depth for the /z"*" 
amounts to 344 /xm with a straggle of 23 ^m. For the 
same energy a proton penetration depth is 74 /im with a 
straggle of 2.3 //m. In this energy range, the deceleration 
mechanism occurs through electronic interactions. As ex- 
pected, the heaviest particles have a shorter penetration 
depth. Figure n illustrates this behavior for particles of 
various masses namely: /x"*", "'^H"'", ^H"^, '^H''", r"*", tt^ as 
possible particles for sample irradiation and ^-^11+, ^'^11+ 
and for hypothetical irradiation. It is striking to find 
a linear relationship on the log-log plot of the penetration 
vs. mass. Additional modeling effort was focused to the 



comparison between penetration depths of /z"*", and 
■^H"'' at various energies. Figure ^ (lower panel) shows a 
comparison between the three particles as a function of 
their energy. 

/iSR Results 

For all samples, transverse-field (TF) /iSR studies were 
first performed in an applied magnetic field of 7 G. This 
low field value was choosen in order to differentiate the 
diamagnetic muon fraction from the muonium one (i.e. 
representing the formation of muonium in the sample). 
For all samples, a reduced fraction (typically 25% of the 
total signal) of diamagnetic muons was observed at all 
investigated temperatures. The spin of these muons pre- 
cesses according to the gyromagnetic ratio of the free 
muon [7^/(27r) — 13.55 kHzG~^ - the muons in a muo- 
nium state exhibit a precession which is typically two or- 
ders of magnitude faster] . The remaining fraction of the 
/iSR signal was related to the occurrence of muonium. 
However, a muonium signal could only be detected at 
very low temperature (below 6-7 K) where it is charac- 
terized by an extremely fast depolarization rate pointing 
to muonium delayed formation. Figure [3 exhibits a typi- 
cal /iSR spectra recorded at 4 K on monoclinic zirconia. 
Whereas diamagnetic muons are detected over the en- 
tire time window, neutral muonium is detected only at 
very early times (bottom panel of Fig. The popu- 
lation of prompt muonium is, however, rather low and 
reaches about 10% at this temperature, A slight tem- 
perature increase has a devastating effect on the prompt 
muonium population and already at 6 K the population 
basically disappears (< 1%), most probably in reason of 
delayed muonium formation It is consequently impossi- 
ble already at low temperature to study the diffusion of 
muonium, even in a semi-quantitative way. Therefore, 
diffusion studies have been restricted to the diamagnetic 
muon fraction. For such TF /iSR measurements with an 
external magnetic field of 50 G were performed. This 
method was preferred to the zero-field /iSR technique 
since it allows the determination of the temperature evo- 
lution of very weak muon depolarization rates quite ac- 
curately. 

For all samples, the diamagnetic part of the /zSR sig- 
nal could satisfactory be fitted utilising the so-called 
Abragam formula, which describes the muon diffusion 
in transverse field (see for example Ref. [l^): 

G{t) = exp[-M2T| exp(-i/rc)-l+t/rc] cos(wt+0) (2) 

where M2 is connected to the second moment of the field 
distribution at the muon stopping site (M2 = 7^(S^)) 
inside the lattice and is the hopping time (which is 
related to the hopping rate F = 1/tc). The oscillatory 
term, on the right-hand side of Eq. ||2Jl reffects the muon- 
spin precession with w = 'y^Bi (where Bi is the internal 
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FIG. 3: Arrhenius plots of the parameter Tc measured differ- 
ent monoclinic samples of Zr02: in the powder sample (solid 
symbols), in the low-temperature sintered nanoparticle sam- 
ple (pellet sample - open squared) and in a natural sample 
(open disks). Note the increased errors bars for the nanopar- 
ticle sample due to the unavailability of the MORE setup for 
those measurements. 



field at the muon site and of the order of the external field 
of 50 G) and (/> is a phase essentially given by the geom- 
etry of the instrument. The Abragam function describes 
the time-dependence, introduced by the muon diffusion, 
of the interaction between the nuclear spins and the muon 
spin. If the muon motion is fast enough, the muon spin 
will sense an average nuclear field distribution, that is, 
the observed muon depolarization rate will decrease. In 
NMR the similar effect is known as motional narrowing. 
The Abragam function reduces, as expected, to the static 
depolarization function (Gaussian) in the limit of slow 
muon motion. On the other hand, for fast motion, the 
limit is correctly given by the exponential depolarization 
function and an overall decrease depolarization rate. 

In the Eq. Q, the parameter M2 was experimentally 
determined at low temperature, i.e. below about 20 K. 
In this temperature region, the depolarization of the dia- 
magnetic /iSR signal is found to be temperature indepen- 
dent and is interpreted as a regime with a slow hopping 
rate (i.e. F — 1/tc <C l/r^ ), i.e. with a muon essentially 
static. On a second step M2 was kept constant during 
all the remaining fitting procedure, corresponding to the 



assumption that the muon diffuses to equivalent sites. 
The parameter M2 was found to equal 0.0256 MHz^, 

0. 040 MHz^ and 0.23 MHz^ for the pellet, powder and 
natural sample, respectively. Such values are unexpect- 
edly high in comparison to calculated field distribution at 
the possible muon site. Assuming the monoclinic struc- 
ture and a perfect crystal, the field distribution at the 
muon site is produced by the nuclear magnetic moments, 

1. e. principally by the ion ^^Zr^"*". This ion has a spin 
1 = 5/2 with a moment /i = —1.298 //n and a relatively 
low abundance of 11.2%. For the most likely interstitial 
sites, and in particular for the site (0.5,0.5,0.5), theoret- 
ical calculations provide a parameter -M2caic of the order 
0.005 MHz^. This huge discrepancy suggests that this 
specific material has inherently a large density of defects, 
impurities and/or vacancies. Consequently, the samples 
may contain important local distortions, not reflected in 
theoretical calculations. The major role played by de- 
fects on the magnetic field distribution at the muon site 
is illustrated by the rather high value of M2 observed in 
the natural sample. In this sample, it is believed that 
a large contribution to the enhanced AI2 is provided by 
iron-rich clusters as FeO or Fe203. 

Figure 13 exhibits the parameter Tc obtained on the 
powder, pellet and natural samples by fitting Eq. ^ to 
the data. Clearly, the temperature evolution of the Tc 
assumes an Arrhenius behaviour of the form: 

— = i/oCxp(-£;o/fcB7') , (3) 

where Eq is an activation energy and fee is the Boltz- 
mann constant. Best fits of Eq. ^ on the data provides 
activation energies of 3.6 ± 0.3 meV, 4.7 ± 0.2 meV and 
16.2 ± 0.9 meV respectively for the powder, the natural 
and the pellet sample (see also table Pi . Eventhough the 
errors bars are rather higher, note that the hopping time 
Tc extracted from the data obtained in the nanoparti- 
cle sample is rather longer than that of the other sam- 
ples, corresponding therefore to a reduced hopping rate 
r. Such observation may be related to the presence of 
lattice strains leading to a reduction of the hopping prob- 
ability through tunnelling. 

From the value of Tc, and taking into account the 
nearly cubic symmetry of the zirconia system, the dif- 
fusion parameter D can in principle be extracted by the 
expression 

(4) 

Tc 6 

Equation assumes a perfect crystal, where the muon 
diffusion occurs between adjacent and identical intersti- 
tial stopping sites. Considering the most probable muon 
stopping site (0.5,0.5,0.5) mentioned above, the jump dis- 
tance d can be shown to be d = where a represents 
the lattice parameter. However, and as observed above. 
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TABLE I; Comparison of the difltusion coefficient (estimated 
at 500 K) for various species in zirconia. The activation en- 
ergy and species radius are also reported. 
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the picture of a clean and perfect crystal has to be ques- 
tioned in view of the quite strong depolarization rate of 
the muon at low temperature (i.e. when the muon is es- 
sentially static). An alternative view is to consider that 
after thermalisation the muon comes at rest at a trapping 
site, which provides a deeper potential than a regular in- 
terstitial site. Such a scenario appears conceivable if the 
concentration of trapping sites is very high and/or if the 
'free' diffusion is extremely fast, meaning that all muons 
would be likely to have reached a trap during the first 
few nanoseconds. Equation (O would still be valid, but 
now Tc would represent the steps in diffusion between the 
trapping sites and M2 would be connected to the second 
moment of the field distribution at a trapping site. As 
observed, such field distribution can be safely expected 
to strongly deviates (i.e. be wider) than the calculated 
one for normal interstitial sites in a perfect crystal. In 
this 'dirty limit', and by assuming that the concentration 
of trapping sites is less than one site per unit cell, the de- 
termination of the diffusion parameter [see Eq. |0J] can 
be considered as a lower limit. 

Table HI presents a comparison of the diffusion coef- 
ficient D obtained for various species in zirconia. The 
values are estimated around 500 K using, if necessary, an 
Arrhenius law. The diffusion coefficients of ions such as 
and Cs"*" in zirconia are reported in the literature 
[l3|. At room temperature typical diffusion coefficients 
are found in the range 10^^^ to 10~^^ m^s~^ for and 
Cs"*", respectively. The diffusion coefficient of hydrogen 
in zirconia (10~^^ m^s~^) was measured by Aufore |l4| . 
This study suggests a significant mobility in zirconia. On 
the other hand, Lim et al. studied by electrochem- 
istry the mobility of proton in zirconia scales and found 
larger diffusion coefficient than for hydrogen. As already 
discussed, from the present study no information can be 
extracted for the muonium diffusion. However, an evalu- 
ation of the muon diffusion coefficient can be performed 
with the restrictions reported above concerning trapping 
effects. For sake of comparison, the coefficient Tc was ex- 
trapolated at 500 K using Eq. Q and subsequently the 
lower bound of the diffusion coefficient D was estimated 
using Eq. 
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FIG. 4: Comparison between the temperature evolution of 
the depolarization rate measured in different Zr02 samples. 
For sake of comparison, the data were here analysed with an 
exponential decay of the polarization (see also text). A nat- 
ural sample (open circles), an ultra- pure powder (solid sym- 
bols) and a layer sample obtained from a corroded Zircaloy 
plate (open triangles). (Note that the error bars are smaller 
than the data points and that the layer sample data were 
performed during 4 different cooling and warming procedures 
with perfect reproducibility). 



Finally, preliminary /xSR studies have been performed 
on a monoclinic sample present as zirconia layer from a 
corroded Zircaloy sample. In view of the reduced thick- 
ness of the layer (100 /zm) and of the available muon 
energy (about 4 MeV), a series of test runs have been 
performed to define the optimum degrader (460 /xm) . For 
the degrader, ultra-pure Al was utilized. The obtained 
thickness corresponds closely to the one calculated from 
the muon range (see above). As the temperature depen- 
dence of the /iSR signal is strongly non-monotonic com- 
pared to the former samples, an analysis based on the 
Abragam formula [Eq. (0)] is not possible and a simple 
exponential depolarization (reflecting satisfactorily the 
overall temperature dependence) was assumed. Figure 01 
shows the T-dependence of the exponential depolariza- 
tion rate of the corrosion layer sample compared to the 
former samples. The different maxima observed point to 
a number of trapping-detrapping effects occurring in the 
layer. Clearly, this suggests a variety of defects in the 
corrosion layer material. This variety impeeds a more 
quantitative analysis of the data obtained. Hence, each 
trapping site type i can be characterized by a 'mean time 
of stay' Tq and the average time needed for a muon escap- 
ing from this kind of trap to reach the next trap variety 
T^. The depopulation and repopulation processes lead 
to the characteristic humps observed in the temperature 
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evolution of the muon depolarization rate. To compli- 
cate even more the picture, each type of trap has a given 
concentration q and a given local field distribution Mj. 
By adding also the fact that the diffusion out of traps 
processes assume their own Arrhenius law (i.e. with dif- 
ferent activation energy for the Tq), it is obvious that 
quantitative information is, at least, extremely complex 
to extract. The present measurements should be basi- 
cally taken as a strong hint for different trapping site 
types in the layer sample. 

CONCLUSION 

The first goal of this muon spectroscopy study was to 
provide more insight on the hydrogen and proton dif- 
fusion in zirconium oxide. The work was achieved by 
analysing muon behavior in zirconia samples. The muon 
interactions in monoclinic polycrystalline zirconia sam- 
ples were investigated, with the aim of applying later 
the study on thin corrosion layers onto Zircaloy. Prompt 
muonium was not identified in zirconia above 6 K. The 
diffusion coefficient of muon in zirconia was estimated to 
be about 10~^^ m^s^^ at around 500 K. In a prospec- 
tive way the work has been completed by studying the 
properties of muon in a Zircaloy corrosion layer. 
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